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The program objectives are fully defined in the original proposal en- 
titled “A Research Program in Active Control/Aeroelasticity in the JIAFS at 
the NASA Lan 7 ley Research Center" dated August 1, 1961. 

Research Conducted during this report period is described below: 
Development of Synthesis Methodology for Multi-functional Active Control 
System . 

INTRODUCTION 

In the last progress report, formulation of return difference matrix 
singular values and their gradients with respect to control law parameters, at 
the plant output was reported. A software for numerical computation was 
developed. This and previously developed Linear Quadratic Gaussian (LQG) 
analysis software are presently being used to generate an optimal design 
scheme. The design objective is to minimize the system root mean square (RMS) 
response and control activity and also enhance the system stability robustness 
by increasing the minimum singular value of the return difference matrix at 
input and/or output. 

Since the system response and stability robustness are contradictory in 
nature, the constrained optimization approach is expected to arrive at an 
optimal compromise solution. 

Optimization Scheme 

The system state-space equations and block diagram of the system are 

shown in Figure 1 . The external inputs u , w and v may be interpreted as 

com com 

white noise processes in an LQG analysis. 

Performance Index (PI) 

The performance index to be minimized is defined as the expected value 

IT 1 IT 1 

J « Ely Q t y + u Q 2 u ] t+0c (1) 

at steady state condition. The output y^ may be replaced by the design output 
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y D if desired. Here y 1 is used to interpret the PI J in terms of singular 
values to be described later. 

Response Constraints (RC) 

The response constraints are defined as expected steady state value of 
the design responses due to stochastic input (mean square responses). 


» A ® — "-r- - 1 

fv 1 

l "d j 

<0 RMS 


i=1 « » » * Np( 2 ) 


Equation (2) may include RMS sensor outputs, control surface deflection and 
rate, modal deflection and rate, structural bending moment, shearforce, 
torque, etc. The design response constraints and their gradients with respect 
to the controller design parameters are computed by solving a set of Lyapunov 
equations. 

Stability Margin Constraints (SMC) 

The stability margin constraints at the plant input is defined as 

A N 2 

g . , = E [Max{0 , [ a - o(I + KG(s))]}f /N (3) 

n +1 . , -D 

D 

<0 

The ® D is the desired minimum singluar value of the return difference matrix 

at plant input (see position (1) in figure 1) and is directly related to the 

simultaneous gain and phase margins defined in Publication 1. See Publication 

2 for detailed explanation of equation (3). 

The stability margin constraint at plant output is defined similarly as 

A N 

g„ _ - E [Max{0 , [a - a(I + GK(s))]}f /N (4) 

D 1=1 

<0 3=jw^ 

The constraints (3) and (4) and their gradients with respect to the controller 
design parameters are obtained over appropriate frequency range. 

The optimization starts at a stable design point where any or all of the 
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constraints may be violated. The designer can choose any set of constraints 
from equations (2), (3) and (4), to be satisfied. The program computes the 
gradients of PI and the active and violated constraints. The usable-feasible 
direction method of CONMIN optimization program is employed to arrive at a 
feasible design point. 

The optimization scheme is presently being evaluated for a drone lateral 
attitude control system and will later be applied to a large space structure. 
Relation between Performance Index and Singular Values 

To aid choice of performance index weighting matrices Q 1 and Q 2 and the 
noise intensity matrices R^, R^, and R^ of the white noise processes u com 
(fictitious), w and v, the performance index J in equation (1) is expressed in 
terms of singular values of various transfer function matrices and is shown in 
Figure 2. 

To understand the significance of each of the six terms, the general 

input-output transfer-function matrix relations are shown in Figure 3. 

From equations (5) and (7) it is easy to see that if Q =R =R =[0], 

1 v w 

Q 2 =R u =I then minimizing J would improve stability margin at plant input. On 

the other hand, if Q =R =R =[0], Q =R =1 then minimizing J would improve 

2 u w 1 v 

stability margin at plant output. Each term in equation (5) can be 
interpreted by examining the input-output transfer function relations in 
equation (6). 

CONCISE STATEMENT OF RESEARCH ACCOMPLISHED 

A control law synthesis methodology for multifunctional active control 
system to satisfy RMS load and response constraints as well as to meet stabil- 
ity robustness requirements at plant input and output has been developed. 
Modem control theory, singular value analysis and optimization techniques 
have been utilized. All stability and response derivative expressions were 
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derived analytically for sensitivity study. The software is incorporated as 
an update to the AB/LAD general control design software package PADLOCS (see 
publication no. 3). 
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FIGURE -1 BLOCK DIAGRAM 
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general input- output relations 
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